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SURVEY OF LITERATURE ON FLOW CHARACTERISTICS OF A SINGLE 
TURBULENT JET IMPINGING ON A FLAT PLATE 
by James W. Gauntner, John N. B. Livingood, and Peter Hrycak* 

Lewis Research Center 

SUMMARY 

The results of an extensive survey of the flow characteristics of a single jet imping- 
ing on a flat surface are presented herein. The flow field consists of four distinct re- 
gions, each of which is discussed in detail. Methods for predicting the spread of the jet, 
and velocity and pressure profiles in the various regions, are presented. 

Comparisons between theory and experimental data are presented for velocity dis- 
tributions in the various flow regions. Experimental data for velocity fluctuations due to 
turbulence are also shown. Some suggested velocities on which to base heat-transfer 
correlations for use in impingement cooling design and some areas requiring further 
investigations are also discussed. 


INTRODUCTION 

A promising method for cooling turbine vanes or blades is to impinge cool air on 
the internal surfaces of the vanes or blades. Although numerous investigators have 
studied the impingement of air jets on flat plates, heat-transfer results differ so widely 
that no adequate heat-transfer correlation for design use is available. It is for this 
reason that the NASA Lewis Research Center initiated an impingement cooling program. 
The ultimate goals of this program are to develop design heat-transfer correlations to 
apply to cylindrically shaped leading edges of vanes or blades and for impingement cool- 
ing along the vane or blade suction and pressure surfaces. The initial phase is limited 
to impingement on a flat plate. Such results may be applicable for cooling blade suction 
and pressure surfaces. 

Development of adequate heat-transfer correlations depends on a knowledge of flow 
characteristics. A number of investigators have studied the flow phenomena of free jets 
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issuing from both circular holes and slots. Some have also investigated jets impinging 
on flat surfaces. Our purpose herein is to summarize many of these previous studies 
and to include derivations of useful equations for determining pressures and velocities 
in the flow field. This report is based, in part, on a survey of literature on flow and 
heat transfer caused by impingement on a flat plate completed under NASA Grant 
NGR 31-009-004 by Professor Peter Hrycak of the Newark College of Engineering. This 
report is restricted to the study of the flow field of a single jet impinging on a smooth 
flat plate. 


DESCRIPTION OF FLOW REGIONS 

The various flow regions formed by a jet impinging on a solid surface have been 
discussed by previous investigators (e.g. , refs. 1 and 2). In the case of an impinging 
jet, four distinct flow regions can be characterized. These are shown in figure 1 
(symbols defined in the appendix) and are described as follows: 

(1) Region I is the region of flow establishment. It extends from the nozzle exit to 
the apex of the potential core. The so-called potential core is the central portion of the 
flow in which the velocity remains constant and equal to the velocity at the nozzle exit. 



Figure 1. - Characteristic regions in impinging jet flow. 
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(2) Region II is a region of established flow in the direction of the jet beyond the 
apex of the potential core; it is characterized by a dissipation of the centerline jet veloc- 
ity and by a spreading of the jet in the transverse direction. 

(3) Region in is that region in which the jet is deflected from the axial direction. 

(4) Region IV is known as the wall jet region, where the directed flow increases in 
thickness as the boundary layer builds up along the solid surface. 


GENERAL FEATURES OF JETS 

Of fundamental importance in the study of impinging jets (or any jet) is the distinc- 
tion between laminar and turbulent jets. Reference 3 states that the critical Reynolds 
number Re Q (based on nozzle diameter and nozzle exit velocity) which distinguishes 
laminar jets from turbulent jets is about 1000. Reference 4 reports four characteristic 
jet patterns for free jets, namely, 

(1) The dissipated laminar jet, Re Q < 300. In this case, the viscous forces are large 
compared to the inertial forces, and the jet diffuses rapidly into the surrounding fluid. 

(2) Fully laminar jets, 300 < Re Q < 1000. In this case, there is no noticeable dif- 
fusion of the jet into the surrounding fluid. 

(3) A transition or semiturbulent jet, 1000 < Re Q < 3000. 

(4) A fully turbulent jet, Re Q > 3000. 

These results are confirmed in reference 5, which also states that jets are turbulent 
for Re Q > 3000. Reference 6, however, refers to a turbulent jet as having Re Q greater 
than 14 000. The reason for this difference in Re Q at which turbulence begins could not 
be resolved from information contained in the references. The fact that differences have 
been obtained indicates that further work is required to determine the factors affecting 
this characteristic. From a heat -transfer standpoint, the determination of Re Q at which 
a jet becomes turbulent is vital because applications of cooling with impinging jets in- 
volve Re Q as low as 3000. 

With regard to different regions of the jet, reference 7 states that two different 
modes of heat transfer exist. The first mode applies to the potential core region, where 
turbulent fluctuations cause the variation in the heat -transfer coefficient. The second 
mode is for the region of established flow, where both turbulence and velocity cause the 
variation. 

The preceding discussion relating jet pattern to Re Q range dealt only with jets 
issuing from circular openings. However, it appears to apply also to jets issuing from 
slots. Reference 8 indicates that slot jets are turbulent for Re Q = 2000 (here Re Q is 
based on the slot width, whereas for circular jets it is based on jet diameter). 
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The disturbance created by a solid wall on an impinging jet is propagated upstream 
at a rate equal to the difference between the sonic velocity and the fluid velocity. The 
wall effect would therefore diminish with increasing Mach number and would finally dis- 
appear when the Mach number equaled 1 and the disturbance created by the wall would be 
unable to propagate upstream. 

In a theoretical study reported by Levey (ref. 9) the effect of the impingement sur- 
face on the jet was not felt by the jet as close as only 1 slot width away from the surface; 
this study was for an idealized case (incompressible flow) with no account for viscous 
effects or turbulent mixing. An experimental study reported in reference 10 confirms 
the theoretical results; no surface effect on the jet was observed beyond about 2 jet 
diameters from the plate. Consequently, an impinging jet can be considered to behave 
like a free jet except in the immediate vicinity of the impingement surface. 


REGION OF FLOW ESTABLISHMENT (REGION I) 

Near the nozzle opening, a mixing zone originates and turbulent exchange of heat and 
mass takes place. The mixing zone grows in width in the downstream direction of the 
jet, leaving a potential core wherein fluid properties and velocity are relatively constant. 
Eventually, the core is dissipated as a result of the mixing action. 

According to Schlichting (ref. 11), by a judicious selection of the origin in the coor- 
dinate system, the half width b (see fig. 1) of the free jet is directly proportional to the 
distance along the jet centerline; that is, 

b = Constant x (1) 

This equation may be used as an approximation for the underdeveloped jet. According 
to reference 12, the pressure in region I is virtually constant and equal to the pressure 
in the surrounding fluid. Hence, the total momentum per unit time should be the same 
for all cross sections normal to the jet; that is, 

J = j U dM = p j U 2 dA = Constant (2) 

J y o 
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Velocity Profile 


Experimental data of reference 13 suggest that the sole force producing the decelera- 
tion of the jet and the acceleration of the surrounding fluid is the tangential shear within 
the mixing region. An approximate analysis of the mean velocity distribution was pre- 
sented in reference 13 and is summarized in this section (see fig. 2). 



The volumetric flow rate Q past successive normal sections may be written as 



U dA 


(3) 


Since Q varies with longitudinal distance x, its ratio to Q q (at nozzle exit) is 


/ OO 

U dA 

J^ = _U 

Q 0 " U o A o " 


(4) 


where A Q is nozzle exit area. Also, since momentum flux J may be expressed as the 
integral of the volume rate of flow U dA times the longitudinal component of momentum 
per unit volume pU, 
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U 2 dA 


( 5 ) 





Also, energy flux equals the integral of the local volume flux U dA times the kinetic 

9 

energy per unit volume p(U /2), 




«o A o 


( 6 ) 


Since the tangential shear within the mixing zone causes the deceleration of the jet 
and the acceleration of the surrounding fluid and since this is an internal process, the 
momentum flux must be constant at every normal section; that is, 


/ 


U 2 dA 


= 1 


U o A o 


(7) 


If viscous action has no influence on the mixing process, the mean flow should be dynami- 
cally similar under all conditions. Thus, a similar velocity profile must characterize 
every section in the diffusion region. As a matter of fact, experimental data follow the 
general trend of the Gaussian distribution 


U 


U. 


= exp 


m 



( 8 ) 


where y is the transverse distance from the jet centerline or axis and a is the standard 
or root-mean-square deviation (the value of y for which U = 0. 606 U m ; see fig. 3). 

This equation is merely a result of a curve fit and has no theoretical basis. The condi- 
tion of dynamic similarity requires that at all cross sections, regardless of the efflux 
velocity, 



( 9 ) 
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for the slot jet and 

-=K 2 (10) 

x * 

for the circular jet (i. e. , the angle of jet diffusion is a constant). Values of K ^ and Kg 
must be determined experimentally. In real jets, the angle of jet diffusion is not neces- 
sarily a constant immediately downstream from the nozzle. In reference 13, a/x proves 
to be different for both the developed jet (region II) and the developing jet (region I). Still 
an attempt has been made in reference 13 to treat the jet as if the same experimental con- 
stant were applicable for both regions. 

Using equations (8) and (10) in evaluating equation (7) yields an expression relating 
the average nozzle exit velocity to the axial velocity at any point downstream of the po- 
tential core for the circular jet 


i 

I 



n 2 tr2 2 

W n 


= i 


u o A o 


( 11 ) 


This equation can then be evaluated at the point where the potential core region ends 
(i. e. , where x equals x c and U m still equals U Q ). The appropriate value of A Q 
may also be substituted into equation (11). Similar substitutions can be made for the 
slot jet, where the momentum is put on a per unit length basis. These substitutions yield 


x c _ 1 

B o 


(12) 
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for the two dimensional jet and 


x 


c 


D 


1 

2K 2 


(13) 


for the circular jet. The nominal inner boundary of the diffusion region satisfies the 
equation 


or 




-\ 


> 


(14) 



2K 2 x 

D 


J 


where B is the width of the potential core at a distance x from the nozzle. The velocity 
distribution in the diffusion zone of region I, obtained by combining equations (8) to (10) 
and (14), becomes 



exp 



2(K x x) 2 


(15) 


for a slot of half width B q /2 and 



exp 



+ K 2 X 



2(K 2 x) 2 


(16) 


for a circular jet. Experimental data are compared in figures 10 and 16 of reference 13 
with results calculated from equations (15) and (16), respectively. Satisfactory agree- 
ment between equation (15) and experiment is shown in figure 10 of reference 13. Experi- 
mental data shown in figure 16 of reference 13 deviate considerably from equation (16). 
This is due to the fact that reference 13 attempted to use a single value of ct/x for both 
the region of flow establishment and the region of established flow. A better comparison 
would result by using the value of ct/x associated with region I. 
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Potential Core Length 


It was previously stated that region I extends from the nozzle to the apex of the poten- 
tial core. The potential core length considered in this report is sometimes called the 
nominal potential core length and its determination requires a knowledge of the centerline 
velocity decay in region II. For that reason, some material which in reality belongs in 
region II is presented under region I and then not repeated in region II. The actual poten- 
tial core length, which corresponds to the spacing between the virtual origin of the jet and 
the point where the jet centerline velocity begins to decrease, is not referred to herein. 

If equation (11) is solved for the ratio of the centerline velocity to the average nozzle 
velocity for the circular jet, and a similar equation is solved for the slot jet, the follow- 
ing expressions result: 

U m. y^q 

u 0 f 7 

W 

for the slot jet and 

1 

^m_^2_ 

U o £ 

D 

for the circular jet. These equations maybe reexpressed, using equations (12) and (13), 

as 


(17) 


(18) 



(17) 
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(18) 


U 


m 


c 

D 


U o * * 

D D 


The nondimensional potential core lengths and Cg in equations (17) and (18) are for 
the slot jet and the circular jet, respectively. 

Although equations (17) and (18) are only valid beyond the potential core (region I), 
they may be used to experimentally determine the length of the potential core in the fol- 
lowing way. Log -log plots of the measured centerline velocity ratio U m /U Q against 
x/B q (or x/D) yield lines with slopes -1/2 or -1, respectively, for the slot jet and the 
circular jet. Extension of these lines to intersect U /U = 1 then gives the desired 
length of the potential core, Cj for the slot jet or Cg for the circular jet. 

Schlichting (ref. 11) presents the following equation for the velocity distribution in 
the axial direction of a developed laminar circular jet, from which an analytical expres- 
sion for the nominal potential core length can be obtained: 


u = li 

8 77 ux 


('*H 


(19) 


where the kinematic momentum flux 


K = 2?r 


f 


U 2 y dy 


and 


=1 /s 
y i6?i 


3K _y_ 
167T ux 


(20) 


Along the jet axis, £ = 0 and the velocity distribution along the axis reduces to 


u 

m 877 ^ 


( 21 ) 


and 


10 


( 22 ) 


K = 




Therefore, 


U m = 


1 ' U o d2 


87 T l/X 


or 





(23) 


This last equation is similar to equation (18) and shows that the potential core length for 
laminar flow (the coefficient of (x/D)”\ which in this case is (3/32)Re Q ) is a function of the 
nozzle Reynolds number. 

For turbulent flow, the axial velocity is represented by equation (19) if, according 
to reference 11 the laminar flow kinematic viscosity v is replaced by the virtual kine- 
matic viscosity e for turbulent flow, where e has been evaluated as 

e = 0.02 56 Uy 1/2 (24) 

It is shown that yocx and U °c x * (for similar flow). Hence, it is possible to express 

e as 

e = i'yRe 0 (25) 


where y is a proportionality factor. The resulting equation for the jet centerline axial 
velocity for the turbulent flow case then becomes 


or 


U m = 


3 U o p2 
32 yvRe o x 


( 26 ) 


_ 3 /x\ 
U 0 ' 32y \D/ 
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Hence, for turbulent flow, the length of the potential core is (3/32)y and is independent of 
the nozzle Reynolds number as long as flow remains similar and the definition of e holds. 

Another expression for potential core length of a turbulent circular jet may be ob- 
tained from the results of Tollmien. He analyzed the flow in the downstream part of a 
turbulent submerged axially symmetric jet. His solution may be used to get the expres- 
sion shown in reference 12 


U m _ 0. 48 D (27) 

U ax 
o 

This result was .obtained with the use of a similarity variable 


<p = X (28) 

ax 

and numerical integration. From table I in reference 14, when U = U m /2 and y = y^yg* 
(p - 1. 23. From experiments described in reference 12, 

y l/2 = 097 x (29) 


Hence, from equations (27) to (29), 

_ 0. 48D<? ^ g i D ^ 

U o *1/2 ’ x 

Comparison of this result with equation (18) shows Cg - 6. 1 for a circular jet. 

Unfortunately, values of core lengths differ for investigations carried out by dif- 
ferent experimenters. The amount of divergence which different observers report 
generally exceeds the expected experimental error and is more than likely explained by 
the effect of the scale and the intensity of turbulence within the jet. Also the actual 
velocity distribution at the nozzle seems to be a factor. The following two paragraphs 
give the results of some of the previous investigators and some information concerning 
their experiments. 

For slot jets, Albertson (ref. 13) states that the potential core extends 5. 2 slot 
widths. He investigated air flowing through slots of widths 1/4, 1/16, and 1/32 inch 
(0 . 6 3 5, 0. 159, and 0. 079 cm) for a range of Reynolds numbers 
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1. 5X10 3 < Re Q < 7. 2xl0 4 

Schauer and Eustis (ref. 15) investigated a 1/2-inch (1. 27 -cm) slot; air impinged on a 

plate located from 10 to 40 slot widths from the nozzle. The values of Re ranged from 
4 4 O' 

2x10 to 5. 3x10 . The core length of about 5. 5 slot widths used in the analysis presented 

in reference 15 was obtained by averaging two similar values obtained from the data of 

van der Hegge Zijnen (ref. 16). Reference 16 considered slot widths of about 0. 197 and 

0. 394 inch (0. 5 and 1. 0 cm) and a value of Re Q of about 1. 3xl0 4 . Reference 17 reports 

results for various -sized rectangular, square, elliptical, and triangular slots, for 
3 4 

3. 6x10 < Re Q < 8. 8x10 and arrived at an average value of potential core’ length of about 
7 slot widths. 

For circular air jets, reference 13 reported a value for C 2 of 6. 2 after testing 
nozzles 1/4, 1/2, and 1 inch (0.635, 1.27, and 2. 54 cm) in diameter for values of Re Q 
from 2. 2xl0 4 to 8. 7X10 4 . Pai (ref. 18) reported a value for C 2 of 6. 5. From the re- 
sults of Tani and Komatsu (ref. 10), the corresponding value for Cg is about 4. 8; this 
result was obtained from data for air flowing through nozzles with diameters of about 
2. 05 and 3. 27 inches (5. 2 and 8. 3 cm) and impinging on a plate 2 to 15 diameters from 
the nozzle. Poreh and Cermak (ref. 1) experimented with water flowing through nozzles 
3/4, 1, and 1^ inches (1. 905, 2. 54, and 3. 81 cm) in diameter impinging on a plate 
located from 8 to 32 nozzle diameters from the nozzle; they report a value for C 9 of 7.7. 
From data shown on figure 2 of reference 1, the local jet centerline velocity can be ob- 
served to be about 10 percent higher than the average nozzle exit velocity. If the poten- 
tial core length of reference 1 were to be reexpressed in terms of the maximum nozzle 
exit velocity, the numerical value would be more comparable to values obtained by other 
investigators. 

In an appendix to reference 13, Baines examined the results of Albertson more 

closely and concluded that the Reynolds number has an effect on the potential core length. 

4 5 

He recommended values for C 9 of 5 for Re = 1.4x10 and 7 for Re = 10 . 

u o o 

Abramovich (ref. 12), however, asserts that the value of Cg does not vary with Rey- 
nolds number but does vary with the velocity profile. However, since the velocity pro- 
file in turbulent flow may be considered to be a weak function of the Reynolds number, to 
some extent C 9 could vary with Re as Baines suggested. However, theory based on 
Prandtl’s second hypothesis shows that, for developed turbulent flow, the potential core 
length is independent of the Reynolds number (eq. (26)). 

To summarize, it has been noted that, for turbulent flow, potential core lengths vary- 
ing from about 4. 7 to 7. 7 have been obtained by different investigators. . Moreover, some 
investigators state that core length is dependent on Re Q , while others state that core 
length is independent of Re Q . This latter view is supported by equation (26) of this re- 
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port. Since, from a heat -transfer standpoint core length is important, it appears that 
additional work in this area is required. 


Centerline Turbulence Distribution 

Investigations have shown (refs. 6 and 7) that the local heat-transfer coefficient at 
the stagnation point on the impingement plate reaches a maximum value at the apex of 
the potential core. Since the centerline velocity is constant in the potential core, the 
observed increase in heat -transfer coefficient must be due to the effects of turbulence. 
Reference 19 shows that the turbulence level on the jet centerline begins to rise at an 
axial distance of only 1 diameter from the nozzle mouth; this is shown in figure 4, taken 
from reference 19, for both a 1-inch (2. 54 -cm) and a 3 -inch (7. 62 -cm) heated circular 
jet. Nozzle size was also shown to affect turbulence and heat transfer. Consequently, 
any reference velocities that may be used or applied in heat-transfer analysis or correla- 
tion may need to be modified to include the effects of turbulence. 



Figure 4. - Turbulence distribution along jet centerline. 
(Data from ref. 19.) 
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Pressure Distribution 


The study of turbulent free jets by use of the boundary -layer theory assumes negli- 
gible pressure variation in the flow direction. However, reference 20 reports obser- 
vation of a radial pressure gradient at a value of x/B Q = 2 (within the potential core) 
with the maximum pressure at the centerline higher than ambient. A minimum pres- 
sure lower than ambient was observed away from the centerline of the jet. Reference 20 
also reports a static pressure ridge, higher than ambient and decaying rapidly with x, in 
the potential core. The ridge is attributed to the persistence of the pressure distribution 
created within the nozzle. Outside the core in region I, all measured static pressures 
were below the ambient pressure. 

From available theory and data, it appears that the effect of a pressure gradient in 
the potential core zone is relatively unimportant with respect to flow development, but it 
seems to warrant an additional experimental investigation. 


Secondary Effects 

The preceding discussions were limited to a study of the length of the potential core 
and the variation of velocity and pressure in the diffusion region. Secondary effects 
produced by the impingement surface may also have some effect on the jet. These 
secondary effects are (1) a dynamic instability effect, and (2) a blockage or interference 
effect. The former is characterized by a lateral oscillation of the complete jet flow 
field while the latter is distinguished by changes in core length and jet spreading rate. 


REGION OF ESTABLISHED FLOW (REGION II) 

Some of the analyses that were required to determine the potential core length and 
the velocity profiles for region I are only applicable to region n. These analyses are 
not repeated, but reference to the appropriate equations is made. 


Centerline Velocity Decay 

Based on discussion contained in reference 13, the centerline velocity decay for the 
slot jet is given by equation (17) 
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(17) 



and for the circular jet by equation (18) 



(18) 


Many experimenters (e. g. , Albertson et al. (ref. 13)) have investigated the centerline 
velocity decay, and their results have verified the general form of these equations. 


Velocity Profile 

According to Schlichting (ref. 11), problems in free turbulent flow are of a boundary - 
layer nature and can be studied with the aid of equations of boundary -layer type. For the 
two-dimensional, steady -state case, these equations are 


uiU + v 3U = 1 0T 

(31) 

3x 3y p 3y 


l U + 9V=0 

(32) 

ax ay 



where t denotes the turbulent shearing stress and where the pressure term has been 
dropped (to a first approximation, the pressure is assumed to be constant). 

To obtain the velocity distribution from these equations, it is necessary to express 
r in terms of parameters of the main flow. This may be done in three ways: 

(1) Use of Prandtl’s mixing length theory, where 

,2 lau 

r = -puv = pi — 

|9y 

(2) Extension of the mixing length theory, where 


au 

8y 


(33) 
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( 34 ) 


and i and i ^ are to be regarded as purely local functions 
(3) PrandtPs second hypothesis 


r = pe — = pxi b (U. 
ay 1 


-U . )9U 
max min' - 

ay 


(35) 


where b is the half width of the mixing zone (from centerline to outer jet 

boundary) and Xi is an empirical constant. From this, e = x,b (U - U . ) 

ji x max min 

is the virtual kinematic viscosity, assumed to be constant over the whole width 
and previously used in equation (24). 

Tollmien (ref. 14) obtained the velocity distribution by solving the system of equa- 
tions (31) to (33); the use of equation (33) involves the application of PrandtPs mixing 
length theory for evaluating the shear stress. 

Gortler (ref. 21) solved the system of equations (31), (32), and (35); this solution 
involved the use of PrandtPs second hypothesis and, according to reference 11, results 
in a simpler solution. Reference 11 further states that best agreement with experiment 
is obtained when Gortler’ s analysis is considered. 

Tollmien* s solution requires the use of successive approximations. For this reason, 
and because of the preceding statement, this solution is not reproduced herein. Refer 
to reference 14 for details of the solution for the two-dimensional and circular jets. 

In reference 14, values of U/U m are tabulated as a function of y/(xlp), where 

is evaluated as 0. 063 and c is the scaling factor in PrandtPs mixing length theory. 

Instead of using y/(x as the abscissa, it is common practice to plot U/U m against 

y/yi/ 2 , where Y-^/2 * s va ^ ue y f° r which the velocity is equal to U m /2. Hence, 

when U/U = 0. 5, y/y« /0 = 1. For the circular jet, a plot of Tollmien’ s tabulated data 
m x/ u > 

shows U/U m = 0. 5 when y/(x 3 P) = 1. 23; hence, if each abscissa is divided by 1. 23, 
values of y/y-^ are obtained, and the ordinates can then be plotted against these new 
abscissas. The same procedure was used to replot the slot jet solution. This is the way 
Tollmien’ s solutions are plotted in this report. 

For the two-dimensional jet, according to Schlichting, Gcrtler’s solution may be 
written as 


17 



u =y-y^( 1 - tanh2 ’h) 

— 2t ?i " ^nh 2 ~ tanh 77 ^ 


where the constant kinematic momentum flux K = J/p and 


x= 2 


where s is the fixed distance from the nozzle, U is the centerline velocity at s, 

’ m, s j > 

e s = WVs 


V 1= — 


For this case, X is the empirical constant. Its value was obtained as 7. 67 by Reichardt 
(ref. 22 ). 

For the circular jet, Gortler’s solution is 


u= A K 

877 ex 


4 T 7T 


(■•} 


x H* 


where K is kinematic momentum flux (constant), e is virtual kinematic viscosity 
(constant), and . 
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For this case, Reichardt determined the empirical constant to be 62. 1. Equa- 

tion (39) is identical with equation (19) with v replaced by e. 

Corrsin (ref. 19) reports experimental transverse velocity profiles through region I 

for a 3 -inch (7. 62 -cm) diameter heated jet and a 1-inch (2. 54 -cm) diameter heated jet. 

Distributions were given for values of x/D of 2, 4, and 6 for the 3-inch (7. 62 -cm) jet 

and a value of x/D of 5 for the 1-inch (2. 54 -cm) jet. In region I, the velocity in the 

potential core remains constant at U„ = U and falls off in the transverse direction 

mo 

exterior to the potential core (see fig. 2). In order to compare these experimental 
velocity distributions with those of Tollmien and Gortler, it was necessary to translate 
the jet centerline to the border of the potential core at the respective values of x/D, 
determine the value of y^yg (where U = U m /2) for each of the distributions, and then 
calculate the abscissa y/y^ * or se l ec ted values of U/U m for each of the distribu- 
tions considered. Figure 5 shows the comparison of the experimental distributions of 
reference 19 with the Tollmien and Gortler distributions. The good agreement in fig- 
ure 5 indicates that the theories of Tollmien and Gortler are satisfactory for use in pre- 
dicting the velocity profiles in region I of the jet even though these theories were devel- 
oped for region n. 





(b) Distance from nozzle in axial 
direction, x = 4D ; nozzle diameter, D = 3 
inches (7.62 cm). 



Dimensionless transverse distance, y/y^ 


(c) Distance from nozzle in axial direction, (d) Distance from nozzle in axial direction, x = 

x « 6D; nozzle diameter, D = 3 inches 5D ; nozzle diameter, D = 1 inch (2. 54 cm). 

(7. 62 cm). 


Figure 5. - Velocity profiles in circular turbulent jets. 


19 



Under the assumption used in the analyses previously discussed, the transverse 
velocity distributions differ only by a scale factor. At all sections beyond the tip of the 
potential core, the velocity distribution is continuous across the line of symmetry, and 
the ratio a/x is therefore indicative of the rate of outward spread of the jet. By the 
same procedure as that explained for the region of flow establishment (region I), refer- 
ence 13 shows that 


U - " V)7K 1 X \ 2K 2 x 2 


for the two-dimensional jet, and 


_U_ 


* — exp ( — 

2K*x 2 


2K 2 x 


(41) 


(42) 


for the circular jet. Equations (41) and (42) show good agreement with experimental 
data for developed turbulent flow. Equations (41) and (42) were obtained by assuming 
the same experimental constant a/x applied to both regions I and II. In addition, they 
were obtained by use of equation (8), a curve fit of experimental data, rather than from 
the governing equations considered by Tollmien and Gortler. 

Figure 6(a), taken from reference 11, compares experimental data compiled by 
Forthmann (ref. 23) with the Tollmien and Gortler solutions for a two-dimensional jet 
and region II. Figure 6(b), also taken from reference 11, compares experimental data 
compiled by Reichardt (ref. 22) with the Tollmien and Gortler solutions for a circular jet 
and region II. These curves are plotted with U/U m as a function of y/y^/2’ con_ 
version of Tollmien' s solution to this type of plot was discussed previously. The experi- 
mental data agree better with Gortler's solution over most of region n, except farther 
away from the center of the jet, for y/yj/2 " > 

Reference 1 also presents a plot of experimental data together with the theoretical 
curves derived by Gortler, namely, 



(43) 
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Distance from 
nozzle exit 
(ref. 23), 
x, 

in. (cm) 

© 7.87 (20) 

3 13.79 (35) 

© 19. 70 (50) 

© 23.61 (60) 

• 29. 55 ( 75) 

Tollmien (ref. 14) 

Gortler (ref. 21) 


Distance from 
nozzle exit 
(ref. 22), 
x, 

in. (cm) 

• 7.87 (20) 

© 10.24 ( 26) 

O 17.72 (45) 

Tollmien (ref. 14) 

Gortler (ref. 21) 


(b) Circular turbulent jet. 

Figure 6. - Velocity distribution in two-dimensional and circular turbulent jets. (From ref. 11). 


and an error curve 



exp 



(44) 


for the transverse velocity distribution in the region of established flow for a circular jet. 
The measured transverse velocity distribution is well approximated by either Gortler' s 
solution (eq. (43)) or the error curve (eq. (44)). 

In van der Hegge Zijnen (ref. 16), slot jet data are compared against the error curve 
of equation (44) except that the constant is 70. 7 instead of 100. Good agreement is 
achieved. 
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Turbulence Distribution Along Centerline 


Corrsin (ref. 19) measured the turbulent velocities U*/U for a 1-inch (2. 54-cm) 
diameter heated jet at various distances from the nozzle. A plot of these values along the 
centerline in region II resulted in a straight line, on log -log paper, with slope equal to 
-0. 88 and intercept at 0. 84; that is, the equation of the line may be written as 



(45) 


An experimentally determined value of U'/U at a value of x/D = 5 (region I) obtained 
in reference 19 does not fall on the extension of this line; hence, it appears that in re- 
gion I the relation of equation (45) does not apply. 

Figure 7 shows a normalized turbulence distribution along the centerline for a cir- 
cular and a slot jet obtained from references 19 and 20, respectively. Reference 19 
shows a plot of U m /U Q against x/D. By taking the data points for U*/U Q at any 
x/D and dividing by the value of U m /U 0 at the same value of x/D, a curve of U*/U m 
against x/D can be obtained. This plot is shown as the solid curve in figure 7. Also 
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shown in figure 7 are several experimental values of U’/U for a 1/2-inch (1. 27 -cm) 
slot jet (ref. 20). 

In summary, if the length of the potential core is known, equations (18) and (45) can 
be combined to give the relation between U’/U and x/D in region II. 

Pressure Variation Along Jet Centerline 

In the analysis of the free jet, ordinarily no variation of static pressure is assumed. 
By neglecting the turbulent y stress, Tollmien (ref. 14) found the centerline static pres- 
sure to be expressible as 


2_ xt 2 


Pm-ex = 0 - 263I 9> u ; 


m 


(46) 


where Kg is an experimental constant expressing the rate of spread of the jet, where 
the half width of the jet is defined as (u/u 

m ) 2 dy. 

According to Reichardt (ref. 22), Kg = 0. 075. Because of the value of Kg, Tollmien 
concluded that the change in static pressure could be neglected. 

Pai (ref. 18) stated that Tollmien’ s approximation, which assumes no variation in 
static pressure, could be improved. He obtained 


m 


- P, 


0. 00165p 



for the slot jet, and 


P 


m 


"PA 


-0. 0029 5p 



(47) 


(48) 


for the circular jet. These equations indicate an overpressure in the middle of the jet 
for the slot jet and an underpressure for the circular jet. However, Townsend (ref. 24) 
studied the magnitudes of the various terms in the equation solved by Tollmien and con- 
cluded that the lateral gradients of the turbulent y stress and static pressure are the 
dominant terms and that the pressure relation should read 

p + pV ,2 =p x (x) (49) 
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where p^(x) is a function of x equal to the static pressure outside the flow. Equa- 
tion (49) applies for both regions I and n, whereas Tollmien's solution applies only for 
region II. 

The result of Townsend's theoretical study is in direct contradiction to Tollmien's 
prediction for the centerline static pressure for a slot jet (i. e. , eq. (47)). Townsend’s 
result yields an underpressure in the middle of the slot jet. This underpressure has 
been verified by data of Miller and Comings (ref. 20), which show a dependence of x/B Q . 
Downstream of x/B Q = 10, these data can be represented to within 4 percent by the fol- 
lowing equation: . 



0. 016 


(50) 


Reference 16 states the velocity on the jet axis U m can be obtained from the total 
pressure head h m in the jet axis from the relation 


U 


m 




(51) 


where f j is a correction factor for the effect of turbulence. The factor f ^ has been 
calculated on the assumption that the total pressure head is mainly affected by the axial 
component of the instantaneous velocity (U + u) , where u is the velocity fluctuation in 
the axial direction, so that 


f 


1 “ 



(52) 


U' = V u^. 


This result has not yet been confirmed by 


where 
authors know. 

The peak pressure that existed along the jet centerline in 
approaches ambient in region n. 


experiment, as far as the 
the core region eventually 
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Static Pressure Profiles 


Static pressure profiles were measured for a slot jet (ref. 20) and are shown in 
figure 8. The trend in the shape of the curves is toward a stable functional form, vary- 
ing only in amplitude as the distance from the nozzle increases (see eq. (50)). Refer- 
ence 20 describes the two lower curves by the relations 


'P ‘ p a\ 


pu: 


m, 


x/B o =30 


= -0. 054 exp [-0. 3 

b 2 


'P ‘ 


p u: 


m 


x/B o =40 


= -0.060 exp (-0. 3 — 
b 2 


(53) 


(54) 



Ratio of distance from jet axis 
to jet width, y/2b 


Figure 8. - Static pressure profiles in 
fully developed region for a slot jet. 
(From ref. 20.) 


Figure 7 of reference 20 shows a plot of the turbulent x -stress profiles in the fully 
turbulent region. Miller and Comings state that the numerical centerline values on 
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this figure are nearly equal to the empirical constants in equations (53) and (54). These 
equations indicate that a static pressure field which is lower than ambient exists through- 
out the turbulent region of the two-dimensional free jet. 


Spread of Free Jet 

The spread of the free jet can be determined as soon as the velocity profile as a 
function of x and y in region II is known. Different experimenters have investigated 
the spread of the jet. For example, Forthmann (ref. 23) examined the spread of the jet 
where y = y, /9 and obtained a rate of spread of U/U equal to 0. 097. 

If the dimensionless velocity profile behaves as exp (-70. 7 y vx* 5 ), as shown in 
reference 16, the rate of spreading of U/U m = 0. 5 will be 0. 0979, compared with 
Forthmann' s value of 0. 097. Reference 12 presents experimental data of Forthmann and 
others, which show the spreading of a slot jet to be the same as the spreading of a circu- 
lar jet. 


REGION OF DEFLECTION (REGION III) 

Velocity Distribution 

As the flow approaches the impingement plate, flowing in the axial direction, the 
velocity decelerates rapidly and the pressure increases. At the stagnation point, the 
velocity is zero and the pressure is a maximum. If the stagnation point is assumed to 
be the origin, the flow is then in the negative z direction (see fig. 1). For three- 
dimensional, axisymmetric frictionless flow, Schlichting (ref. 11) states that the axial 
component of the velocity may be written 

U = -2az (55) 

and the radial component of the velocity as 

V = ar (56) 

where a is a constant. When z is replaced by z n - x, equation (55) becomes 
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(57) 


U = -2a(z n - x) 


or 


— = - 2/1 -JL\ (58) 

aZ n \ Z J 


This is the form of the equation found in reference 1 for the axial velocity distribution in 
the deflection region; reference 1 found that 


a = 98 \E 


(59) 


where K is the kinematic momentum flux. Similarly, reference 1 found 


V 

az 

n 


_r 

z 


n 


or, since a 


= 98\/k, 


V 

Vk z „ 



(60) 


(61) 


The preceding analysis was obtained for irrotational flow, which, according to ref- 
erence 10, precludes immediate application of the theory to the impingement region. A 
theory was developed in reference 10 for axisymmetric inviscid rotational flow. A stream 
function which automatically satisfied the continuity equation was inserted into the mo- 
mentum equation (with the pressure term eliminated). Axial and radial velocities were 
found to be 


— = -f + -r 2 (f" - b 9 ) - ir 4 K,f + -f””J + . . . (62) 

aL 4 z 8 \ 4 8 / 


_V_ 

aL 



I r 2 f*M 
8 



K,f’ + if*"" 
4 8 



(63) 


where L is the jet radius at a fixed distance from the nozzle, -f represents the axial 
velocity along the jet centerline, and primes denote differentiation witti respect to z; the 
jet centerline velocity is approximated by 
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(64) 


f(z) 


bjZ 



+ — b„z 3 + 
3 , 0 



b 6 z 
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where b^ to bg and are constants. From conditions applied in the analysis, 
bj = 1 and b^ = 0. Values of bg, bg, bg, and bg are evaluated so that the polynomial 
f(z) agrees with measured centerline velocities for small values of z, where inviscid 
flow holds. The constant is determined such that the value of U observed at some 
distance from the plate is expressed by the first three terms of equation (62). The radial 
velocity along the plate can then be expressed as 



(65) 


It is shown in reference 10 that the surface velocity computed by equation (65) shows 
good agreement with the experimental data. 


Pressure Distribution 

Denoting p’ as the total pressure on the jet centerline and p as the static pres- 
sure at any arbitrary point, reference 11 states that 


p’ - p =£ (V 2 + U 2 ) = p — (r 2 + 4z*) 

2 


( 66 ) 


Because this equation was derived for frictionless or irrotational flow, the value of p' 

is constant along the impingement plate and equals p* Equation (66) may be rewritten as 

s 


— - -2- = — (r 2 + 4z 2 ) 
pg Pg 2 g 


(67) 


At the plate, z = 0, so the equation now becomes 


a 2 z 2 / \2 

h ’ - h = (—} 

S 2 S vnj 


( 68 ) 


Reference 1 found that a = 98^K with K = M/p. Therefore, equation (68) may now be 
written as 
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(69) 


From plots of the pressure and velocity distributions in reference 1, it was concluded 
that the stagnation zone could be defined to be interior to a hemisphere with the center at 
the stagnation point and with a radius of 0. 05 z n (i. e. , r/z n < 0. 05 and z/z n < 0. 05). 

The effect of diffusion is important beyond these limits. On the boundary of this hemi- 
sphere, the pressure intensity deviates from the parabolic shape resulting from equa- 
tion (69) and approaches the ambient static pressure with increasing r/z n . The region 
of established radial flow starts where the pressure gradient is nearly zero (i. e. , 
r/z n > 0. 3 (ref. 1)). 

Experimental data obtained for a circular jet are presented in reference 10. The 
pressure difference along the impingement plate (PiO^q - was made dimensionless 
by use of the nozzle dynamic pressure pU^/2. The pressure fell from a maximum value 
p’ at the stagnation point to p. (ambient) at a radius of 1. 6 to 3 nozzle diameters; hence, 
the region of impingement, which was quite small, was found to extend from 1. 6 to 
2. 2 diameters axially and 1. 6 to 3 diameters radially, with the stagnation point at the 
center. 

Velocity profiles were obtained at various radial stations and indicated that maxi- 
mum values existed near the surface. Plots of these maximum velocities and of those 
computed from the measured pressure distribution (assuming no dissipation in the flow 
(i. e. , p^ - (Pr) z _Q = \ ^m ) s h°wed that the measured and computed values coincided up 
to 0. 6 to 1. 2 diameters for nozzle to plate spacings of 4 to 12 diameters, respectively. 
Hence, it appears that the effect of viscous dissipation is negligible up to these distances, 
which, according to reference 10, are about 40 percent of the radial extent of the deflec- 
tion region. 


WALL JET REGION (REGION IV) 

Velocity Decay Along Plate 

In analyzing the so-called wall jet region (region IV, fig. 1), Glauert (ref. 25) divided 
the flow region into two parts: (1) an inner layer where the effect of the wall is present, 
and (2) an outer layer which is characterized by the features of a free turbulent flow. The 
velocity at the boundary between these two layers is the maximum velocity. It has been 
shown (refs. 26 to 28) that this maximum velocity may be expressed by. the relation 
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Constant 


(70) 


V 


m 


n 


r 


Four different values of n (1. 14, 1. 1, 1. 143, and 1. 12) have been found by each of four 
different investigators (refs. 25 to 28). Glauert (ref. 25) was the first, however, to 
formulate equation (70) using n = 1. 14, which gives results very close to experimental 
values. In equation (70) the power of r is greater than 1 for a radial wall jet, whereas 
for a free jet it is equal to 1; hence, the decay of the maximum velocity along the plate 
is faster than that of a free jet. 


Velocity Profile Through Wall Jet 

Glauert obtained solutions for the two layers of flow mentioned previously. For the 
inner layer, he assumed the 1/7 power velocity profile; for the outer layer, he assumed 
that the eddy viscosity was constant. The solutions were then matched at the boundary 
between the two layers, where the shear stress was assumed to be zero. The result was 
a velocity distribution throughout the jet, dependent on a single constant a. Since a 
varies with Reynolds number, the velocity profiles vary with distance. However, the 
variation in a is very small, and so the velocity field is approximately similar. 

Other investigators (refs. 1 and 12) found it useful to employ a reference boundary 
velocity V RB where 


_ Constant 
V RB ~ 


(71) 


The reference velocity V RB is the centerline velocity at r of a hypothetical free jet 
issuing radially from the stagnation point and can be obtained by extrapolating the veloc- 
ity profile of the outer jet region to the wall (see fig. 1). GOrtler’s free-jet solution for 
the outer layer was discussed earlier and can be used as follows: 


V 


V 


RB 


1 - tanh 2 



(72) 


where is an integration constant. 

Abramovich (ref. 12) used the relation 


\ 
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^U(i-C 3/2 ) 2 

V m 


to represent the velocity distribution through the outer layer of the wall jet; here, 


C = 


z - 5 
b’ 


(74) 


Figure 9 shows a plot of V/V m against z / z 1 /2 * or theories of Glauert, 
Abramovich, and Gortler. To obtain the curve for Gortler' s distribution as shown on 
figure 9, equation (72) was multiplied by V^g/V m . K K 5 = 10 (ref. 1) and r = 
(ZjygK-s/Kg), then ^RB^m’ anc * ^6 can 1,6 ca -t cu l a ted by using the boundary conditions 
V/V m = 1/2 at z / z i /2 = t and V/V m = 1 at z/ = 0.2 (ref. 25). The three curves 
in figure 9 show fairly good agreement. 



Figure 9. - Comparison of wall -jet theoretical velocity pro- 
files. 
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Reference 1 reports experimental data for the velocity distribution through the wall 
jet. If a value of = 10 is used in equation (72), good agreement with the data is ob- 
tained. This is shown in figure 10. Here, in the outer zone, the experimental data can 
be approximated almost equally well by 


— = 1 - tanh 2 10 - 


z > 6 


V. 


RB 


or by the error curve 



(75) 


(76) 


Bakke (ref. 28) showed by means of experimental data that the velocity gradient near 
the wall is much greater than that of Glauert’s analysis (ref. 25). Figure 11 shows a 


Slot width, 



Figure 10. - Velocity distribution in region of radial flow (region IV). (From ref. 1. ) 
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Figure 11. - Comparison of velocity profiles of wall jet by 
Glauert and Bakke. 


comparison between Bakke' s data and Glauert' s analysis; it can be seen that, except for 
the profile in the boundary layer near the plate, the similarity law is well satisfied. 

As experimentally determined by Bradshaw and Gee (ref. 29), Schwarz and Cosart 
(ref. 30), and Kruka and Eskinazi (ref. 31) the turbulent shear stress was zero not 
at z = 6, but rather at a point closer to the wall. Also, several investigators (refs. 26, 
29, and 32 to 34) have reported that the turbulence effects were found to be considerably 
higher than those corresponding to a flat plate. These turbulence effects probably contri- 
buted to producing a velocity profile near the wall whose power slope was measured as 
1/11 to 1/14. 

Spalding (ref. 35) found a method where no knowledge of the shear stress at z = 6 
is required, and McGahan (ref. 36) showed that this method and the two -layer method 
with shear stress zero at z = 5 give essentially similar results. 

More recently, Cadek (ref. 37) analyzed the behavior of two-dimensional impinging 
jets, based on the assumptions of normal impingement, a flat, constant temperature 
impingement surface, incompressible flow and constant fluid properties, and wall 
temperature higher than jet temperature. 
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The impingement region treatment includes Cadek's own static pressure measure- 
ments; integral methods are applied to the calculation of the inner flow associated with 
the wall -boundary layer. The outer layer is assumed to behave as a free jet. 

The analysis is based on integral methods appropriate for the laminar flow near the 
s tag nation region. An approach using the boundary -layer thickness at the position where 
the velocity within the boundary layer has been reduced to just 50 percent of the maximum 
velocity in the boundary layer has allowed Cadek to eliminate the effect of the distance 
from the jet orifice x/B Q . 

In summary, the integral approach using the two -layer model, despite its inconsis- 
tencies, provides a generally good means of predicting velocity profiles. Additional 
investigations are in order to obtain more information about the power slope of the pro- 
file near the wall-, assumed as 1/7 for theory and obtained as 1/11 to 1/14 for data. 


Spread of Wall Jet 

The spreading of the wall jet in terms of the half value can predicted by 

use of the equation 


-1I2l = Constant ^ (77) 

The values of the exponent m were found in references 26 to 28 to be 0. 9, 1. 028, and 
0. 94, respectively. 

The spread of the wall jet, as in the spread of the free jet, can be determined as 
soon as the velocity profile as a function of r and z is known. If the dimensionless 
velocity profile behaves as given in equation (76), the rate of spreading with respect to 
r of the half width of the wall jet (where V/V m = 0. 5) will be 0. 083. Experimental data 
in figure 7 of Poreh et al. (ref. 26) show the rate of spreading of the wall jet to be 
approximately 0. 09. 

POSSIBLE APPLICATIONS FOR USE IN HEAT-TRANSFER CORRELATIONS 

The ultimate objective of the impingement cooling program was the determination of 
heat -transfer correlations adequate for design purposes. Although this report is re- 
stricted to a study of the flow characteristics of an impinging jet, it seems in order to 
make a few remarks relative to what information may be useful in obtaining the desired 
heat -transfer correlations. 
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Since the centerline jet velocity is constant in the potential core and decreases 
beyond the apex of this core, maximum heat transfer will be obtained if the impingement 
plate is separated from the nozzle exit by the length of the potential core. If this length 
is not known, an assumed value of about 6. 1 nozzle diameters may be used. In determin- 
ing the heat transfer at the stagnation point, the use of the arrival velocity in evaluating 
the Reynolds number has proved useful (refs. 6 and 7). The arrival velocity is that 
velocity which would exist at the same distance from the nozzle as the stagnation point 
if no impingement plate were present. 

If the plate lies inside the potential core, the arrival velocity is equal to the nozzle 
exit velocity since the velocity within the potential core is constant. If the impingement 
plate lies outside the potential core, the arrival velocity may be calculated from equa- 
tion (18) 


U C 9 

-2 l=-l (1 

U o £ 

D 

where x is replaced by z r , and U m then becomes the arrival velocity. 

Similarly, the velocity fluctuations due to turbulence may be calculated from equa- 
tion (45) 


tt . /v \-°.88 

— = 0. 84 /— J (45) 

u 0 W 

A combination of equations (45) and (18) permits calculation of U* in terms of 
U , as previously illustrated in figure 7. If heat -transfer data do not correlate on the 
arrival velocity, this velocity may be corrected by use of the fluctuation velocity, and 
another attempt at correlation may be made. 

To supply adequate impingement cooling, the proper spacing of free -jet nozzles must 
be known. This nozzle spacing is dictated by the rate at which the jet spreads. 
Abramovich (ref. 12) reported a value of cKyjygJ/dx = 0. 097. 

Since the primary interest in the impingement flow investigation is the determination 
of heat transfer, some remarks relative to heat transfer at the stagnation point are in 
order. In order to estimate the heat transfer at the stagnation point of an impinging 
flow, the local radial velocity gradient (dV m /dr) r _Q, which appears as a parameter in 
the usual stagnation -point heat-transfer equation, must be evaluated. This can be done 
experimentally by relating this parameter to the static pressure distribution on the sur- 
face in the immediate neighborhood of the stagnation point. Assuming incompressible 
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flow, this parameter is proportional to the square root of the curvature of the pressure 
distribution at the stagnation point and can be obtained from experimental data from the 
equation 



where 

r w wetted radius of impingement surface 
T stagnation point temperature of fluid 
R gas constant 

p ’ total pressure at stagnation point 

This parameter can be calculated as follows: to determine this parameter, pressure 

data in the form 1 - (p/p') are first plotted as a function of (r/r )^. The slope of the 

s w 

resulting curve is determined graphically at r = 0, and the value of (dV m /dr) r _Q is then 
computed from the preceding equation. 

At distances along the plate from the stagnation point, it has been shown that a 
maximum velocity exists at the boundary between an inner flow layer, where wall effects 
are felt, and an outer layer, which is characterized by the features of a free jet. This 
velocity is expressed by equation (70) 



for which values of n ranging from 1. 1 to 1. 143 have been reported. 

If heat -transfer data in the wall -jet region cannot be correlated by use of the maxi- 
mum radial velocity, other radial velocities must be used. Three methods for predict- 
ing the velocity profile through the wall jet have been discussed. 

Although the present review was restricted to single jets impinging on a flat plate, 
the information presented may be applied to the use of multiple jets. Since the spread 
of the jet is discussed, this information may be applicable for locating other jets in an 
array so that maximum heat -transfer coefficients could be achieved along a surface. In 
this way, for example, it may be possible to use rows of jets to obtain maximum cooling 
along the suction and pressure surfaces of turbine vanes or blades. 
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There is some evidence (ref. 38) that heat transfer from impinging jets on flat plates 
and on curved surfaces is similar. With additional flow studies of jets impinging on 
curved surfaces to verify similar flow characteristics, the discussions presented herein 
may eventually be applicable to such surfaces as turbine vane or blade leading edge 
regions. Another investigator who has studied heat transfer from the inner side of con- 
cave surfaces is D. E. Metzger (ref. 39). This study is directly applicable to impinge- 
ment cooling of vane or blade leading edges. 


CONCLUDING REMARKS 

An extensive survey of the flow characteristics of a single jet impinging on a flat 
surface led to the following conclusions: 

1. The flow field can be divided into four distinct regions. Methods for predicting 
velocities and pressures in the various regions agree reasonably well with experimental 
data. 

2. Additional work is required to determine the influence of flow Reynolds number on 
core length and the large discrepancy between the theoretical and experimental velocity 
gradients near the wall. 

3. The length of the potential core was found experimentally to vary from about 4. 7 
to about 7. 7 nozzle diameters by different investigators. This length has also been shown 
theoretically to equal 3/32 times the nozzle exit Reynolds number for a laminar jet and 

to be independent of nozzle exit Reynolds number for a fully developed turbulent jet. By 
combining Tollmien’s analysis for a fully developed turbulent jet with an experimental 
value for the spreading of the half width of the jet, the potential core length was found to 
be about 6. 1 nozzle diameters. 

4. The spreading of a free jet where the velocity equals one -half the maximum veloc- 
ity was found to be 0. 097 times the distance from the jet nozzle. 

5. The axial velocity for a free jet is constant and equal to the nozzle exit velocity 
over a distance equal to the potential core length. 

6. The axial velocity downstream of the potential core for a circular free jet is in- 
versely proportional to the distance from the jet nozzle expressed in nozzle diameters. 
For a slot jet, the axial velocity downstream of the potential core is inversely propor- 
tional to the square root of the distance from the jet nozzle in slot widths. The constants 
of proportionality are the potential core length for the circular jet and the square root of 
the potential core length for the slot jet, respectively. Impinging jets behave like free 
jets up to several nozzle diameters from the surface. 


37 



7. An equation for determining velocity profiles normal to the jet axis through the 
potential core is available; it depends on a single experimental constant, namely, the 
length of the potential core. This equation was developed by Albertson. 

8. Velocity profiles beyond the potential core region can be determined from the 
anasyses of Gortler and Tollmien. Equations for these profiles depend on a single 
empirical constant. Good agreement with experimental data has been achieved. 

9. Equations are available for calculating turbulence velocity fluctuations for the 
circular jet if the length of the potential core is known. An assumed potential core length 
of about 6. 1 nozzle diameters may be used if the actual length is not known. 

10. The radial velocity component in the deflection region is proportional to the 
radial distance from the stagnation point. 

11. Three analytical solutions (those of Abramovich, Glauert, and Gortler) for the 
velocity profiles through the wall jet region are available. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 5, 1969, 

720-03. 
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APPENDIX - SYMBOLS 


A area of jet cross section 

a constant 

B width of potential core 

B q slot width 

b half width of free jet 

b* distance from edge of 

boundary layer to edge 
of wall jet 

bp. . .,bg constants in polynomial 

approximation for veloc- 
ity along jet centerline 

Cj dimensionless potential 

nominal core length 
for slot jet 

Cg dimensionless potential 

nominal core length 
for circular jet 

c constant 


D nozzle diameter 

E energy flux 

f poly nominal approxima- 

tion for velocity along 
jet centerline 

f^ turbulence correction 

factor 


f’,. . . ,f m ” derivatives of f with 
with respect to z 

g acceleration due to 

gravity 

h pressure head, p/pg 


h’ 


J 

K 

K 

L 


l, l x 


M 


m, n 

P 

P’ 

Q 

R 


Re„ 


s 

T, 


total pressure head, 

P’/pg 

total momentum flux 

kinematic momentum 
flux, J/p 

. ,Kg constants 

jet radius at fixed dis- 
tance from nozzle 

mixing lengths 

mass 

exponents 

static pressure 

total pressure 

volumetric flow rate 

gas constant 

Reynolds number based 
on either nozzle diam- 
eter or slot width 

distance from stagnation 
point in radial direction 

fixed distance from nozzle 

stagnation point tempera- 
ture 


t 

U 

U» 


u 


time 


velocity in axial direction 


time mean of square of 
instantaneous velocity 


fluctuation ,Vu2 


instantaneous axial veloc- 
ity fluctuation 
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V 


V’ 


'RB 


v 


X 


y 


* 1/2 


z 


n 


z l/2 


a 

y 

6 

e 

? 


^1 


velocity in radial direc- 
tion 

time mean of instanta- 
neous radial velocity 
fluctuation 

reference velocity 

instantaneous radial 
velocity fluctuation 

distance from nozzle in 
. axial direction 

absolute potential core 
length 

distance from centerline 
in radial direction 

value of y for which 
U = 1/2 U 

m 

distance from impinge- 
ment plate in axial 
direction 

distance between nozzle 
and impingement plate 

value of z for which 

V = 1/2 V m 
m 

constant 

proportionality constant 

boundary -layer thickness 

virtual kinematic vis- 
cosity 
z - 6 
b 

Xy' 

x 



v laminar flow kinematic 

viscosity 


% 

P 

a 

T 

<P 

*1 

Subscripts: 

m 

max 

min 

o 

r 

s 


sc 

w 

X 



16;r ux 
density 

standard or root -mean - 
square deviation 

turbulent shearing stress 


y_ 

ax 

empirical constant 


centerline or maximum 

maximum 

minimum 

average nozzle exit 

radial position 

stagnation point value or 
fixed distance from 
nozzle 

settling chamber 

wetted 

ambient 
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